Spectrum of n Electrons in Graphene as an Alternant Macromolecule 
and Its Specific Features in Quantum Conductance 
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An exact description of it electrons based on the tight-binding model of graphene as an alternant, 
plane macromolecule is presented. The model molecule can contain an arbitrary number of benzene 
rings and has armchair- and zigzag-shaped edges. This suggests an instructive alternative to the 
most commonly used approach, where the reference is made to the honeycomb lattice periodic in 
its A and B sublattices. Several advantages of the macromolecule model are demonstrated. The 
newly derived analytical relations detail our understanding of n electron nature in achiral graphene 
ribbons and carbon tubes and classify these structures as quantum wires. 



I. INTRODUCTION 

The electronic spectrum of 2D graphite, one-atom 
thick hypothetical material with the structure of a hon- 
eycomb lattice, was first described more, than sixty years 
ago in the band theory language [1]. Since then it was 
addressed by many authors, especially after the discovery 
of multi-wall [2, 3] and single- wall [4, 5] carbon tubes and 
free standing, experimentally accessible graphene [6-10]. 
Two approaches have dominated theoretical modeling of 
graphene and its daughter structures [10]. One following 
the line of the Wallace pioneer work produced a con- 
siderable development based on the Solid State theory 
methods [11-13], see also [14, 15] and references therein. 
The other approach has focused on peculiar electronic 
properties near the Fermi energy, where the methodology 
of quantum electrodynamics proved to be both heuristic 
and instrumental [16]. 

The molecular approach to the description of graphene 
has received comparatively little attention. Zigzag car- 
bon tubes can be thought as build up of cyclacenes [17]; 
similarly, a graphene sheet can be constructed from lin- 
ear acenes. From this point of view, graphene is a typical 
alternant macromolecule [18-20] representing a vast field 
which has been contributed in a number of fundamental 
studies [18, 19, 21-24]. 

Recently, simple tight-binding model of n electrons in 
graphene as a macromolecule was briefly reported [25] . It 
was demonstrated that such an approach provides a deep 
insight into the graphene electronic structure. Several 
developments have been outlined which were inaccessi- 
ble or missed in the previous studies. In particular, a 
controversy between the obvious nonequivalence of arm- 
chair and zigzag directions in the honeycomb lattice and 
the identity of band structures near the Fermi energy for 
armchair and zigzag carbon tubes is resolved by exposing 
principal distinctions between different types of metallic 
graphene ribbons and related carbon tubes in accurate 
analytical expressions. 



The purpose of this article is to give a detailed de- 
scription of the model proposed in [25] and partly, its 
further developments [26-29] . In Section II supplemented 
by Appendixes A and B, the solution of the Schrodinger 
problem for the graphene sheet framed by armchair and 
zigzag edges is obtained. It comes to the dispersion rela- 
tion, where wave vector components are not independent 
but interrelated via a transcendent equation, typical for 
alternant oligomers. The usage of this pair of equations is 
specified for graphene daughter structures, armchair and 
zigzag carbon tubes and parent graphene ribbons, as the 
basis for further analysis. The Fermi energy region re- 
ceives much of attention in Sec. Ill and Appendix C, 
where it electron spectra of achiral graphene ribbons and 
carbon tubes are expressed in terms of elementary func- 
tions. As a particular application of these results, quan- 
tum conductance of graphene-based wires is discussed in 
Sec. IV. This is followed by a concluding Section. 



II. EXACT SOLUTION OF THE 
SCHRODINGER PROBLEM 

Shown in Fig. 1 is a honeycomb NxAf lattice. 
The lattice label indicates that in the armchair direc- 
tion, the sheet of graphene contains N hexagons in 
polyparaphenylene-like chains, whereas in the zigzag di- 
rection, it has Af hexagons forming acene chains. Hydro- 
gen atoms along edges are not shown and not taken into 
account in the nearest-neighbor tight-binding Hamilto- 
nian H. By exploiting m,n, and a = l,X,p,r labeling 
explained in Fig. 1, the ir electron wave function that 
satisfies the Schrodinger equation, 

H^ = E^, (1) 

can be represented in the form of expansion 

\m,n, a), (2) 

m n a=l,r,\,p 



* aleonSif m .liu.se 



where |m, n, a) is the 2p z orbital at the ath atom of ben- 
zene ring with coordinates {m, n}, the summation is run- 



2 




FIG. 1: Indication of labels of carbon atoms used in the 
present description of the n electron spectrum of N x J\[ hon- 
eycomb lattice, a is the minimal translation distance in the 
lattice. 



ning over all sites of the honeycomb lattice, 



EM ij . Trim 



JV+1 



, sin 



7rj(m-l/2) 

aT+i 



a = L r, 



a = A, p, 



(3) 



1, 2, A/", are to be found from 



and unknowns 4>> n a , j 
Eq. (1). 

Details are described in Appendixes A and B, where 
solutions to Eq. (1) are obtained for the open boundaries, 
that is for the lattice terminated by armchair and zigzag 
edges as they appear in Fig. 1. For periodic boundary 
conditions in x direction, the wave function is also found. 

In Appendix A, the tt electron spectrum is shown to be 
determined by equation (energy is in units of the hopping 
integral \t\ [1, 14]) 



Ef 2 = l±4cos 



TT] 



2(Af- 



T) COS 2 



3 +4 cos 2 



TTJ 



2{N+l) 



, (4) 



where for the sign + and — , and each value of j, N values 



of Kj , «• Vj = 0, 1, N-l, are solutions to equation 



sin K f(Ar + l/2) 



= =f2cos 



TTJ 



2(Af+l) 



(5) 



Thus, tt electron states of graphene can be classified 
into 2Af "j-minus" and "j-plus" conduction and equal 



structures of graphene, the spectrum is completely de- 
termined by an appropriately modified Eq. (4) or its 
substitute Eq. (A15) in the case of periodic boundary 
conditions in y direction. 

In a graphene strip with N » Af, zigzag edges affect 
a small part of the spectrum. In the limit N — > oo, 
these edges play no role at all and can be replaced 
by a continuous variable, the wave vector in units of the 
inverse periodicity V3a in armchair direction, see Fig. 1 . 

Hence, the spectrum of an infinitely long graphene 
strip with armchair edges, called henceforth armchair 
graphene ribbon (aGR), is fully determined by a single 
equation that reads 



< \/3fc| < TT 



(6) 



Here and in what follows, we are using dimensionless 
wave vectors in units of a -1 . 

As shown in Appendix A, a pair of equations similar 
to Eqs. (4) and (5) can be obtained for periodic bound- 
ary conditions in y direction. From Eq. (A15), it follows 
that the spectrum of zigzag carbon tube with an infinite 
length [or, simply, zigzag carbon tube (zCT)] is deter- 
mined by 



:T (fc,±) = ±^/l±4|cos#| 
0,1,...,JV-1 



cos^i+4cos 2 ^. 



(7) 

In the next case, complications of finding the spec- 
trum connected with the necessity to solve Eq. (5) are 
avoided due to that the boundary conditions at zigzag 
edges are chosen to be periodic. In combination with 
armchair open edges, this corresponds to a segment of 
armchair carbon tube (aCTS), the spectrum of which 
can be written as 



E fj' = E )r 



,aCTS 



±yi±4cos2(^ij jcos^f I +4C0S 2 27 
j = l,2,...,Af,f = 0,1,...,N-1 



KA/+T) ■ 



(8) 



In the limit Af — > oo, one can obtain the band spectrum 
of armchair carbon tube (aCT) by making a replacement 

■KJ 

aT+T 



-^3— ky in the above equation. This yields 



into 2jV "j-minus" and "j-plus" conduction and equal aCT \ fe ± ^ ., fc ± 

number of valence bands with N levels (k^,, ) within E j' ) = ±\/l ± 4cos -|- cos +4cos 2 -|- 



each band. This is ANN of the total number 2N(2Af+ 1) 
of tt electron levels. Additionally, there are two A^-fold 
degenerate levels with energies ±1. These originate from 
the states with zero wave-function amplitudes at the I 
and r sites in linear acenes [28]. 

Equation (5) makes one quantum number dependent 
on the other. It appears because of zigzag-shaped edges. 
Before finding solutions to this equation, let us consider 
related structures, where zigzag edges either do not ex- 
ist or their effect can be disregarded. For such daughter 



< k± < TT 



(9) 



Finally, in the limit N,Af — * oo, when both quantum 
numbers in Eq. (4) can be treated as continuous vari- 
ables, the dispersion relation for tt electrons in the infinite 
graphene sheet reads 

E{k±,k±) = ±\J l±Acos^r cos Acqs 2 ~%. 

< V3fc± , fc± < TT 

(10) 




FIG. 2: Lower part: 2D graphite lattice build up of A and B triangular lattices (on the left) and NxAf honeycomb lattice 
of graphene sheet (on the right). Zero-energy points, six for 2D graphite and two for graphene, are indicated as K and cross 
points, respectively. Circle (0, -^-) and star (ttj, ^-) indicate additional zero-energy points of zCT dispersion, Eq. (7), and 
aCT dispersion, Eq. (9), respectively. Distance in k space between cross points (= ^p) differs from that is between K and K' 
points which is equal to Upper part: dispersion in graphene [on the right, Eq. (10)] and in 2D graphite [on the left, Eq. 
(11)]; mid panel shows the Brillouin zone in 2D graphite. Energy is in units of |t|, wave vector is in units of a -1 . 



Usually, another form of dispersion relation is referred 
identically. Obtained for the honeycomb lattice consist- 
ing of two periodic triangular lattices A and B [1] , it reads 
[14, 15, 30-32] 

E w (k x ,k y ) = ±\[l + 4cos^cos^ +4cos 2 \ . 
< k x ,V3k v < 2tt 

(11) 

The dispersion calculated according these two energy- 
wave vector relations, one for NxAf and the other for 
periodic honeycomb lattice, is shown in Fig. 2. The 
striking difference between the dispersion in the left and 
right panels originates from the different symmetry of 
locally identical lattices. 

The subsequent discussion is primarily focused on the 
energy region near the Fermi energy (for this model, it 
equals zero). In the majority of related studies, this re- 
gion associates with zero-energy points of the Wallace 
dispersion relation (11). Therefore, it is worthwhile to 
compare such points in Eqs. (10) and (11). These are 
shown in Fig. 2 and indicated by crosses and by K, 
K', respectively. K points correspond to the corners of 
the first Brillouin zone of 2D graphite (middle panel in 
Fig. 2). Cross points have different coordinates and the 
number of these zero-energy points is only two. The dis- 
tance between cross points in the k space differs from 
that is between neighboring K and K' points. All this 



is a consequence of the lower symmetry of NxAf hon- 
eycomb lattice in comparison with its counterpart, 2D 
graphite lattice periodic in A and B sublattices. 

For carbon tubes, zero-energy points of dispersion in 
zCT, Eq. (7), and in aCT, Eq. (9), are, respectively, 
cross/circle points and cross/star points. The dispersion 
of aGR (6) has only one zero-energy point, it is the cross 
point, and for the zGR dispersion, this point cannot be 
shown in the real k x -k v plane. Different manifestations of 
zero-energy points, as summarized in Table 1, in the band 
spectra of achiral carbon tubes and graphene ribbons will 
be illustrated below by concrete examples. 



III. BAND STRUCTURE NEAR THE FERMI 
ENERGY 

In this section, we take a closer look at the spectra 
of graphene daughter structures specified above, namely, 
finite- width graphene strips which are infinite in armchair 
(aGR) and zigzag (zGR) directions, and corresponding 
carbon tubes which can be viewed as zGR ruled in arm- 
chair direction (aCT) and aGR ruled in zigzag direction 
(zCT). Note that in this description, periodic boundary 
conditions along the graphene ribbon and respective car- 
bon tube are not used. Such boundary conditions are rel- 
evant to a GR ring and CT toroid, the structures which 
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TABLE I: Coordinates of zero-energy points (z.e.p.) in k 
plane for achiral carbon tubes and graphene ribbons, (a) for 
zigzag graphene ribbons, zero-energy point cannot be defined 
in real coordinates k x and k y ; values of the corresponding 
longitudinal wave vector (in units of a -1 ) at zero-energy point 
(or points) are shown in column ko. gZ is degeneracy of the 
uth band energy \EZ(K = 0)|; (b) (TV + l)/3 is an integer; 
(c), (d) Af is even, odd. ££_ L denotes the HOMO-LUMO gap 
in carbon tube segments (CTS) and graphene ribbon strips 
(GS), see text. 



z.e.p. 

x 



k x , ky ko g v 



2zl 



aCT 



zCT 



2tt 
3 



2tt 

o, 



2tt 
3 



l(i/=0) 
2(^0) 



(b) 

\/3(AA+l) 



2_7T 

3 



(c) 
-2Jv|ln 4 



1(d) 



aGR x 
zGR 



— 



0,^ 



1 2e- 2JV l ln ATTTl 



(a) 




FIG. 3: Band structure of armchair carbon tubes as calculated 
from exact equation (9), solid lines, and its approximations 
(12), dashed lines, and (12) without the q term in parentheses, 
dotted lines; N = 20. Only few conduction bands are shown. 
The valence bands with the same band index are just a mirror 
reflection of conduction bands to negative energies. Spectrum 
refers to the cross and star points; star point itself does not 
belong to the spectrum, see text. 



have not the same spectra as the counterparts with open 
ends [33]. 

For carbon tubes, approximate but sufficiently accu- 
rate expressions will be derived from the exact dispersion 
relations, ft will be demonstrated that in comparison 
with previously suggested approximations, these expres- 
sions improve agreement with the exact model results and 
expose essential differences between the spectra of aCT 
and metallic zCT. For graphene ribbons, the obtained 
results extend and correct earlier descriptions. Also, it 
seems advantageous that the spectra of finite and infinite 
graphene ribbons and carbon tubes are described here on 
equal footings and with the use of the same basic equa- 
tions. In the case of carbon tubes, our approach leads 
to the same spectra as were obtained in the framework 
of zone folding technique [11, 14]. However, the quoted 
technique is not applicable to graphene ribbons that has 
motivated the appearance of other calculation schemes 
to be discussed later on. In comparison with our previ- 
ous related publication [27], this section places a special 
emphasis on the origins and values of band degeneracy 
in graphene ribbons and carbon tubes. 

Below, only a part of the spectrum is in focus. It cor- 
responds to the minus branch of Eq. (4). Therefore, we 
omit indication "— " in energy and wave vector notations. 
Index a will be used for system labeling, a = aCT, zCT, 
and so on. 



series near the cross and star points. Up to the lowest 
order in q = k y — and q, 7r ^ A l7 J - , the minus 

branch of Eq. (9) takes the form 



N 



EZ(q) = ±^-JmZ 2 (l-4q)+q 2 



aCT, 



to,, = 



2ir\v\ 
V3N 



« 1, \q\ «1, v = 0,±1, 



(12) 

Note that Eq. (12) represents expansions near both zero- 
energy points. The state with v, q = is not degener- 
ate because the star point does not belong to the aCT 
spectrum. All other band energies in this spectrum are 
twofold degenerate. Electron states with energies other, 
than the band bottoms or tops, are two- and fourfold 
degenerate for i/=0 and v^O bands, respectively 

Conduction (valence) band bottoms (tops) are at- 
tained at 



V3 



aCT 2 



and are equal to 



aCT 



V3 



,aCT 



-ill 



, 3 . 

1 ra\ 

32 1 



CT 2 



(13) 



(14) 



A. Spectrum of Armchair Carbon Tube 

To obtain an expression of aCT energy spectrum near 
the Fermi energy, Eq. (9) needs to be expanded in power 



The quantity |j/| _1 -^to^; ct would be the band spacing, 
if linear in q term under the root in Eq. (12) were disre- 
garded, as, e.g., in Refs. [15, 30, 31, 34]. 

In the spirit of analogy to be discussed in Sec. 3E, 
we refer to the approximation E°(q) — ml 2 + q 2 
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wave vector ( a 1 ) 



FIG. 4: Full spectrum (on the left) and its view near cross 
point (on the right) for armchair carbon tube (aCT) and 
zigzag graphene ribbon (zGR). N = 20; solid lines are cal- 
culated according Eq. (9); dashed lines represent approxi- 
mation (12) for aCT and approximation (30) for zGR. Only 
conduction bands are shown. 



as relativistic-like approximation. The band structure of 
aCT calculated with the use of exact Eq. (9) and approx- 
imately, according to Eq. (12), and in the relativistic-like 
approximation just mentioned, is represented in Fig. 3. 
Equation (12) provides a very accurate reproduction of 
exact results which need no comments. In contrast, the 
spectrum in the relativistic-like approximation shown by 
dotted lines is noticeably worse. However, both approx- 
imations are practically equivalent for the description of 
quantum conductance, and integral characteristics such 
as density of states. 

According to its derivation, spectrum (12) refers to 
the cross and star points. However, the star point does 
not belong to the spectrum. Thus, with an account to 
spin degeneracy, the electron state with energy E (0) 
is twofold degenerate, and states with energies E^ Jq(0) 
are fourfold degenerate. 

The same conclusion follows from zone folding [14, 15, 
30-32]. However, the band spectrum E v (k y ) shown in 
Fig. 4 associates with K and K' valleys, and not with the 
cross and circle points. In our opinion, this is misleading. 
The distance between these points is different and, as dis- 
cussed above, K points are characteristic for 2D graphite 
and not for carbon tubes and graphene ribbons. 



B. Spectrum of Zigzag Carbon Tube 

An analytical expression for zCT spectrum in the vicin- 
ity of cross and circle points can be obtained from Eq. (7) 
in a way similar to the above consideration. We note 
first that as known, distinct from armchair carbon tubes 
which are always metallic, zCT has a gapless spectrum, 
if j* = Af/3 is an integer, E^,(k x =0) = 0. Other- 
wise, zCT spectrum has a gap. If Af/3 is not an integer, 
band index of the lowest conduction (highest valence) 
band can be equal either to j* = J ^^ 1 or to j* = J ^^ 1 ■ In 
general, these two possibilities correspond to different 
types of semiconducting tubes, although they have the 
same band gap. Distinctions between, let us say N— 1 
and Af+1 semiconducting carbon tubes (and their parent 
graphene ribbons) are seen in the electron states, which 
are obtained in Appendix B, and also, in the respective 
Green's functions, where these distinctions are especially 
prominent [28]. 

Again, Eq. (7) has to be expanded near zero-energy 
points in powers of k x and 27r( ^ J - . In a metallic zCT, 
j* = 4r and j* — 2^ for the cross and circle points, 
respectively. In a semiconducting zCT, the choice of j* 
is explained above. As a result of the expansion, we arrive 
at 



EZ(k x ) = ±&y/mZ 2 + kl 

< = ^( 1 + ^=)« 1 ' - = zCTm, (15) 
ml = \v- §| « 1, cr = zCTs, 

i/ = 0,±l,..., 

where extensions in labeling indicate metallic (zCTm) 
and semiconducting (zCTs) zigzag carbon tubes. 

Distinct from the aCT spectrum is that the v=0 band 
in zigzag carbon tubes is twofold degenerate, as all other 
7T bands near the Fermi energy. Another distinction is 
that within one valley associated with either cross or 
circle zero-energy point, it is only one value of k x that 
matches any energy value within the given band. In the 
aCT spectrum, there are two values of k y which corre- 
spond to the same energy within the uth band. 

It is worthwhile noting that the band energies in the 
zGR spectrum contain a linear in vjM term, whereas 
in the aCT spectrum, the correction term is quadratic 
in v/N, see Eq. (14). Hence, the band spacing in aCT 
spectrum is nearly constant, whereas in zCT spectrum, 
it is essentially dependent on the band index. 

From the comparison of calculations according to Eqs. 
(7) and (15), as illustrated in Fig. 5, it is apparent that 
retaining the v term under parentheses in the definition 
of mf Tm gives a much better agreement with the exact 
calculations than the use of m^, CTm = ^jf^-- The extra 
zz-dependent term in the band energy E^_^ (0) and the 
shift of band bottoms/tops with the increase of band in- 
dex towards larger wave vectors in aCT are missed in 
the k-p approximation which is equivalent to the use of 
(2+l)D quantum electrodynamics formalism [31]. As al- 



6 



zCT 




ft f 

J\r = 


33 

TJ _ _ , 


. p 










-2 




-i. 






+1 






Ji-""" 







0.2 



0.4 



0.6 




FIG. 5: Same calculations as in Fig. 3, but for zigzag carbon tubes, semiconducting TV = 32, 34 and metallic Af — 33. Solid 
and dashed curves labeled by v values represent Eqs. (7) and (15), respectively. Dotted curves in the mid panel correspond to 



zCTm _ 1t\v\ 



Spectra refer to the cross and circle points in Fig. 2. 



ready mentioned, this difference becomes insignificant in 
the limit of large tube diameter. 

An important conclusion that follows from Eqs. (12), 
(14), and (15) is that near the Fermi energy, carbon tubes 
can have three types of spectra: (i) with (to a good 
approximation) equally spaced bands, as for aCT, (ii) 
with regularly irregular band spacing, as for zCTm, and 
(iii) with an alternating band spacing, approximately, be- 
tween values of and as for zCTs. 

Subsequent sections reveal a very close similarity be- 
tween the band spectra of armchair and zigzag carbon 
tubes and their parent graphene ribbons, zGR and aGR, 
respectively. Also, one prominent exception from other- 
wise identical regularities will be described in detail. 



C. Spectrum of Armchair Graphene Ribbon 

The spectra of graphene ribbons have been persistently 
studied numerically at different levels, from tight-binding 
to ab initio [13, 20, 35-40]. Few attempts to attack the 
problem analytically [40-43] have been paralleled by ex- 
tensive simulations of ribbon spectra. Two forthcoming 
sections provide a fully analytical description of aGR and 
zGR spectra near the Fermi energy. This description de- 
tails and extends the results of works [25-27]. 

It is easy to show that with minor changes, which 
are connected with different boundary conditions in the 
transverse direction, aGR spectrum repeats zCT spec- 
trum (15). Specifically, the condition of metallicity for 
aGR requires j* — 2(Af + l)/3 to be an integer. Then, 
index v = j—j* = corresponds to the zero-energy band. 
If 2(7V+ l)/3 is not an integer, aGR spectrum has a gap, 
and the band closest to zero is either j* = (2Af + l)/3 or 
j* = (2Af + 3)/3 depending on which of these two num- 
bers is an integer. There is only one zero-energy point in 
the spectrum of metallic aGR, it is the cross point. 

For k x ,u/M « 1, the exact spectrum (6) simplifies 



to 



< = ^( 1 + i7#TI)) <<1 ' - = *GRm, (16) 
ml = jjZpj- \v - || << 1, a = aGRs, 

i/ = 0,±l,..., 

which repeats the spectra of metallic and semiconduct- 
ing zigzag tubes, where Af 2 (Af +1). Another way of 
derivation of aGR band structure was presented in Ref. 
[42]. It was concluded that "... there is no general rule 
of the subband index The above equation construc- 
tively opposes this statement. 

The difference between the boundary conditions for 
aGR and zCT results in about two-times smaller band 
spacing in the aGR spectrum, than it was found for the 
zCT spectrum. Continuing this comparison we note that 
all v bands in aGR spectrum are nondegenerate. This 
conclusion just repeats the known result [13]. The k-p 
approximation, where electron states in K and K' valleys 
have to be admixed, prescribes the twofold band degen- 
eracy in armchair graphene ribbons [41]. 

D. Spectrum of Zigzag Graphene Ribbon 

Of four basic graphene-based wires, achiral graphene 
ribbons and carbon tubes, only in zigzag ribbons it occurs 
that longitudinal and transverse motions are not separa- 
ble. Analytically, this interrelation was expressed for the 
first time in Ref. [41]. It has the form of a transcendent 
equation 

q = k x cot(V3k x N), (17) 

where q and k x are, respectively, the longitudinal and 
transverse wave vector components referring (according 
to the authors of the quoted paper) to one of K points. 
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Equation (17) was obtained by exploiting the Dirac equa- 
tion for massless fermions with the dispersion 



E{k x ,q)=±^-^kJ 



(18) 



As a particular case, Eq. (17) appears in a recent analy- 
sis of graphene structures with more complex edges than 
just zigzag or armchair [44]. However, the results of these 
studies are restricted by limitations of the long-wave ap- 
proximation. Here, the problem is addressed on the basis 
of exact analogues of Eqs. (17) and (18), represented in 
the next two equations. 

The part of zGR spectrum that includes the Fermi 
energy is described by the minus branch of Eq. (10), 



E (k x , ky) 



±v 1 — 4 COS -j- cos ■ 2 



4 cos 2 



(19) 

where for each value of < k v < ir, allowed values of k x 
has to be found from 



sm V3k x N „ k v 

— = 2 cos -f. 

sin V3k x (N + 1/2) 2 



(20) 



Note that Eq. (18) is just an expansion of radicand (19) 
near the cross point in powers of k x and q = k y — ^ up 
to k x and q 2 . Approximation (17) follows from the exact 
relation (20) after the following replacements 2cos(7r/3 — 
q/2) 1 - V3q/2 and 



sin 



V3k x N 



sinV3fcx(iV+l/2) 



=> 1 



cot 



N 



this cannot be justified in any rigorous way. 

Equation (20) has N solutions k x u \ v = 0, 1, ...,N— 1, 
for each value of k y . One of these solutions is imaginary, 

y/3k x ^ — iS, if k y falls into the interval ^§-+q c < ak y < it, 
as illustrated in Fig. 6. Shown in this figure are graphical 
solutions of Eq. (20) for two cases, q < q c and q > q c . 

q c corresponds to = 0, that 



The critical value of q ■■ 
is 



q c = 2 arccos 



( N 



V 2N + 1 



2n 

T' 



(21) 



In the state with k x 
equals 



and q — q c , the electron energy 



E 



0, 



2tt 



E c 



1 — 2 cos 



(22) 



For energies \E\ < E c , implying that q > q c , the spec- 
trum is described by Eq. (19), where \^3k x = i5, and by 
equation 



sinhcSiV 



sinh5(iV+ 1/2) 



cos ■ 



sm - 



q > q c . (23) 



In this representation too, the quantum numbers q and 
6 are interdependent. The rest of the spectrum is deter- 
mined by Eqs. (19) and (20). 




FIG. 6: Graphical solutions of Eq. (20) which correspond 
to intersections of horizontal lines, representing the right 
hand side, with curves yi (k x ) = sin ^^^ 1/2) an d y 2 (S) = 
siiih^AT+i^) ' N = 50. Blue line intersects only yi; red line 
intersects t/i and yi- 



An equivalent and in certain respects more convenient 
expression of the spectrum can be obtained by combining 
these two equations in order to exclude k y from the first 
of them, see Appendix C. Then, the energy of electron 
states becomes a function of one variable k x . 



E G *{k x ) 



± 



i(V3k x /2) 



iV3k x (N +1/2) 



(24) 



where k x is dependent on k y via Eq. (20). 

Up to this point, our consideration of zGR spectrum 
was exact. Now, we proceed with useful approximate ex- 
pressions which make apparent the electronic structure 
of zigzag ribbons and expose its similarity and dissim- 
ilarity with the band spectrum of armchair tubes. In 
this discussion, the ribbon width will be assumed large, 
N » 1, so that q c = ^ {N \ 1/2) and E c = (2N + 
The approximate formulas to be presented give rather 
accurate estimates already for N ~ 10. 

The dispersion for edge states (ES), for which q > q c 
and \E\ < E c , is described exactly but implicitly by 



E ES (k x =iS/VS) = ± 



sinh(<J/2) 
sinh5(iV + 1/2)' 



(25) 



One can see that for N6 > 1, the edge-state energy goes 
to zero exponentially with the increase of 6, 



E (5) = ±2sinh((5/2)e" 



-S(N+l/2) 



s»i 



■4 ±e 



-SN 



(26) 



The behavior of E (S(q)), that is the edge-state disper- 
sion shown in Fig. 7, is more complicated. Its explicit 
expression can be found from Eq. (23) under certain re- 
strictions on q. 

For 1 < q/q c < 2, this dependence is reasonably de- 
scribed by S w q [l — 2 exp(— 2v3giV)] , so that 



E ES ( q ) 



±V3qey:p{-V3qN). 



(27) 
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FIG. 7: Edge state dispersion in a JV=10 wide zigzag ribbon. 
Solid curves represent exact calculations according Eqs. (20) 
and (25). Various approximations are marked by equation 
numbers. Curve (17), (18) and curve (27) are undistinguish- 
able in semi-log scale. Inset: same data for small values of q 
in non-logarithmic scale. 



This approximation is close to exact calculations up to 
q w n/12. For larger wave vectors, ir/6 < q < n/3 (or 
NS > 2), the solution to Eq. (23) can be represented as 



2 In 2 sin 



tt/3- 



(28) 



showing that with the approaching of the wave vector to 
its maximal value, q — > ir/3 (k y — > it), the edge state 
energy goes to zero as 



E (q) = ±e 



-2iV|ln(f-g)| 



(29) 



Thus, dispersion of edge states, which in numerical 
studies appears as a dispersionless band [13, 36], is gov- 
erned by exponential dependencies defined in Eqs. (27) 
(for q c < q < tt/12) and (28) (for tt/6 < q < tt/3). The 
crossover between the two can approximately be desig- 
nated to a region, where 7r/12 < q < tt/6. Analytical 
expressions of edge state dispersion, which are defined in 
the above equations, are compared with the exact results 
in Fig. 7. As seen, the long-wave approximation followed 
in Refs. [41, 44] and illustrated by curves (17)-(18) and 
(27) does not reproduce the dispersion of edge states in 
the larger part of ES band. 

Out of the ES band, \E\ > E c , and for k x , \q\ « 1, 
Eq. (19) can be rewritten as (see Appendix C [45]) 



EZ(q) = ±^JmZ*(l-&q)+<r i 



i/ = 0, 1, ... << N, 
7^+ 1/2) 

" ~ V3A ' 



cr = zGR, 



where q < q c within the v = band. Without the linear 
in q term, the spectrum (30) can be obtained from Eqs. 
(17) and (18). 

As shown by many numerical calculations (but never 
proved analytically), with an exception of the ES band 
energy interval, zGR and aCT have a very similar band 
spectra. Both spectra have nearly equally spaced bands 
with the bottoms/tops which are pronouncedly shifted 
towards larger wave vectors with the increase of band in- 
dex. The band spacing in zGR spectrum is two times 
smaller, ^37, than it is in the aCT spectrum. There- 
fore, aCT bands are in-between pairs of zGR bands. The 
band structure parameters for zGR, band energies and 
positions of bands are given by the same equations as 

aCT 

for aCT, Eq. (13) and Eq. (14), where m v should be 
replaced by ra v . All these features can be seen in Fig. 
5 that compares the full aCT and zGR spectra of ir elec- 
trons, as well as their details near zero-energy points. 

So, for \E\ > E c , the zGR spectrum can be obtained 
from the aCT spectrum by the "mass" scaling, and the 
other way round. In the narrow energy interval 2E C near 
the Fermi energy, aCT and zGR spectra are qualitatively 
different. This is connected with the appearance of ES 
band in zigzag graphene ribbons. As shown, the dis- 
persion within this band is exponential and goes to zero, 
when the wave vector approaches its maximal value. The 
k-p approximation, as it was used in Ref. [41], leads to a 
"reversed" picture, where ES energy tends to zero when 
the wave vector goes to zero. 

To wind up this section, we note that only marginal 
changes are required to apply the above scheme of cal- 
culations to carbon tubes and graphene ribbons of finite 
length, taking into account their termini. In particular, 
the gap between HOMO (the highest occupied molecu- 
lar orbital) and LUMO (the lowest unoccupied molecular 
orbital) is equal to [46] 



= 2 < 



0, 

2V3(JV+1) : 



e~ 2N \ ln M oddM, 

-2/V|ln -rl 



— ^ integer, 
otherwise, 

even N 1 



a = aCTS, 



a = zCTS, 



= GS, 



(30) 



(31) 

where N and N is the length of aCT and zCT segments 
denoted as aCTS and zCTS, respectively; GS abbreviates 
armchair (N » Af » 1) and zigzag (Af » N >> 1) 
graphene strips. HOMO-LUMO gap values and other 
characteristics which have been discussed here, are sum- 
marized in Table 1. 
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E. Relativistic Analogy 

The band structure of carbon tubes and graphene rib- 
bons can be summarized in a single line [48], 



(32) 



where k° has the meaning of dimensionless wave vector, 
to£ (only here, in units of is the z/th band energy, 

and index a specifies the structure and the ^th band 
minimum/maximum in the ID k space k° as follows 



ir\v-l/5\ 

— 77 — ' 



2n\u-l/3\ 

— ' 

2ir\v\ 

V3N ' 



o, 
o, 

o, 

0, 



,zGR\ 



f+^« CT ) 2 , 



(33) 

a = aGRm, aGRs, zGR, zCTm, zCTs, aCT, from top 
to bottom. Equation (32) is valid for fc£,ra£ << 1; v = 
0, ±1, ... for all structures, except zGR in which case v = 
0,1,... . 

One can see that within each band (except the ES 
band) the dispersion relation (32) has the form of ID rel- 
ativistic energy-momentum relation in its conventional 
representation with the speed of light equal to unity 
[47]. Thus, from the point of view of energy and mo- 
mentum conservation laws, quasiparticles in the u=Q 
conduction/ valence bands of metallic carbon tubes and 
graphene ribbons (except zGR) should behave as ID 
massless Dirac fcrmions (or neutrinos/antineutrinos), 
whereas in bands they should behave as relativis- 
tic particles having mass 

Perfect penetration of electrons/holes with linear dis- 
persion into a classically impenetrable region has been 
noticed in several publications [49-51]. More recently 
[29], it was shown that by passing from the exact tight- 
binding description to approximation (32), expressions 
for the probabilities of tunneling through a potential step 
and through a potential barrier within the ^th band (and 
inter-band scattering prohibited) exactly coincide with 
formulas, derived almost eighty years ago for tunneling 
of relativistic particles [52, 53]. It is really amazing that 
by predicting the penetration probability for massive rel- 
ativistic particles, Oskar Klein predicted the probabil- 
ity of transmission of charge carriers through n/p junc- 
tions in alternant macromolecules and, in particular, in 
graphene, a material unknown in his time. 

Equations (32) and (33) make obvious that spectra of 
achiral carbon tubes and graphene ribbons can be di- 
vided into three groups: (i) metallic spectra with equally 
spaced bands, case of aCT and zGR; (ii) metallic spec- 
tra with an irregular band spacing, where in case of zCT 
(aGR) it can be any fraction of ^ ( a/'+i ) anc ^ al so > larger 
than that; and (iii) semiconducting spectra with a band 



spacing, alternating between Jjjj- and and between 
^ and 3 ff +1 in cases of semiconducting zigzag tubes 
and armchair ribbons, respectively. These differences in 
the band structures along all cis and all trans carbon 
chains reflect in a quantitative manner the nonequiva- 
lence of armchair and zigzag directions in the honeycomb 
lattice structures. In the infinite graphene sheet these 
differences disappear. 



IV. QUANTUM CONDUCTANCE 

Within the framework of the Landauer approach [54- 
57] , the zero-temperature ohmic conductance of an ideal 
wire is equal to 



G%E)=G Y J 9l T uiE), 



(34) 



where Go = 2e 2 /h is conductance quant, g° is the band 
degeneracy (spin degeneracy 2 is included into Go), and 
transmission coefficient T° is zero or unity, depending on 
whether the i^th band is open or closed for charge carriers 
with energy E. 

With an account to Eq. (32), T?(E) = Q(E - ^m£) 

for conduction bands, and T?(E) = Q(\E - ^m£|) for 
valence bands; Q(x) is the Heaviside step function. The 
values of g^ CT and g z J^ represented in Table 1 must be 
doubled because electron/hole states with ±fc^ CT ^ 
and ±k z J$ ^ are degenerate. This "rule" was earlier 
noticed in the numerical study of GR conductance [58] . 

The electron/hole conductance of armchair and zigzag 
carbon tubes and their parent graphene ribbons has thus 
the form of a ladder, symmetrically ascending with the 
increase of energy for electrons, and with the decrease 
of energy for holes. For the charge carrier energy that 
falls in-between the ^th and (z/+l)the bands, the wire 
conductance equals 



G a (E) = Go 



v aGRm, aGRs, 

1v + 1 zGR, 
2v zCTm, zCTs, 

2(2i/+l) aCT, 



(35) 



The width of steps repeats the band spacing m° defined 
in Eq. (33) in units of ^\t\. The height of the z4h ladder 
step is G times band degeneracy as explained above. 
In the aforementioned study [58], the step width was 

calculated numerically. For G" R (E), the same expres- 
sion was suggested, whereas G *" (E) according to the 
quoted paper is equal to 2^Go; for G° (E) no expres- 
sion was presented. The possible reason for the extra 
factor of two in the conductance of metallic armchair 
ribbons is explained bellow. 

As long as dispersion (32) is valid, the conductance 
ladders of basic graphene wires can be classified in the 
following three types: (i) with regular step width for al- 
ways metallic aCT and zGR wires; (ii) with irregular step 
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width for metallic zCT and aGR; and (iii) with alternat- 
ing step width for semiconducting zCT and aGR. These 
characteristics distinguish graphcne ribbons and carbon 
tubes as ID quantum wires from the conductance ladders 
known for 2DEG channels [59]. 

The appearance of G a {E) depends on the energy scale 
which, in its turn, is determined by the ribbon width 
(tube circumference). Experimentally, this can result in 
different observations, from ladders, which might be fully 
or partly resolved, to unresolved ladders. In the latter 
case, that is when the discrete behavior of quantum con- 
ductance cannot be resolved, all ladders are smoothed 
out into straight lines with the same slope and zero- 
energy values equal to G (0) = Go and G (0) = 2Go 
for graphcne ribbons and carbon tubes, respectively. 

Returning to the aGRm conductance, we remind that 

in this case, the step width equals (l+ ^jr^ ) i see 
Eq. (33). Therefore, for large N, quite many pairs of 
bands near the Fermi energy have very close energies, 
i.e., they are apparently degenerate. Then, G* (E) = 
2v, if v « N. Obviously, under the same conditions 
G sCTm (£) = 4za 



V. SYNOPSIS 

A new methodology of analytical modeling of tt 
electron spectrum of graphene and its daughter lattices, 
achiral carbon tubes and graphene ribbons, has been 
presented. It is based on the tight-binding model 
of graphene as a macromolecule with armchair- and 
zigzag-shaped boundaries. The exact solution of the 
Schrodinger problem, the spectrum and wave functions, 
have been obtained and illustrated by a number of ex- 
amples. Several spectral features, which were previously 
accessible only for numerical calculations, have received 
an adequate analytical description in terms of elementary 
functions. In comparison with the understanding based 
on the 2D graphite band structure, the macromolecule 
model is shown to be more relevant and more benefi- 
cial. It provides a consistent description of graphcne 
electronic properties. Presenting the full details of this 
model sheds light on the intimate interrelation between 
graphene, acenes, and other conjugated oligomers. In 
general, this model gives more comprehensive picture of 
what may be called rclativistic appearance of graphcne. 
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APPENDIX A: SPECTRUM 

By neglecting the overlap between tt orbitals at neigh- 
boring C atoms, (m, n, a\m' , nf, a 1 ) — 8 mm i5 nn i5 aa i , Eqs. 
(1) and (2) can easily be transformed into a set of four 
equations 

A + VWl.n.A, (Al) 

(A2) 
(A3) 

+ l/>m+l,n, P - (A4) 

For the NxAf sheet of graphcne, these equations are to 
be solved with open boundary conditions which read 



*Pm,0,r = i>m,N+l,l = 0. 



(A5) 



Coefficients ip m ,n,\ and ip m ,n,p can be expressed in 
terms of ip m ,n,l and ip m ,n,r as 

(E 2 -l)lp m . n ,\=E {lp m ,n,l+4>m-l,n,l)+1pm,n,r + 1pm-l,n,r, 

(A6) 

(E 2 -l)tj) m . n _ p =E {4> m! n,r+4>m-l,n,r) + 4>m,n,l + 1pm-l,n,l- 

(A7) 

Exploiting these two equations in (Al) and (A4), we ob- 
tain a reduced set, involving only ip m .n,i and ip m .n,r- This 
set has the form 

E (E 2 -^)^m,n.l = {E 2 -l)V> m ,n-l,r+-E {tym,n,l+1pm-l,n,l 
+ 1pm+l,n,l) +2lpm,n,r+lpm-l,n,r+1pm+l,n, ri (A8) 

£(£ 2 -l)VW,r = (£ 2 -l)VVn+:U+£ (2VV«,r+^m-l,n,r 
+ VWl,n,r) +'24>m,n,l+1pm-l,n,l + 4>m,+l,n,l- (A9) 

Now, it is convenient to represent wave function coeffi- 
cients in the form (3). By substituting it into Eqs. (A8) 
and (A9) and performing standard algebra, we arrive at 



i,a = S&,i<P n -i,T + 0( n+ i),z, a = l,r, 



(A10) 



where g\ v = = 4 cos\^/2)Vj\ = jfe, g^ = 
g* tr = E[E 2 - 1 - 4cos 2 (^/2)]DT 1 , an d 

Vj = [E 2 -E - 4cos 2 (^/2)][S 2 +E- 4cos 2 (^/2)]. 

(All) 
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kT^J fcT^I 



N 



r I 



M-monomer 

MHM-— -M 



Iff* 
Y3a 



-M 



FIG. 8: NxAf graphene sheet in Fig. 1, represented as a 
sequence of dashed-frame boxes, illustrates the analogy with 
M-oligomer structure exemplified by iV-long oligomer of poly- 
paraphenylene. In Eq. (A10), g\ plays the role of monomer 
Green's function matrix element referring to monomer bind- 
ing atoms I and r. 



The latter equation is nothing else but the determinant 
of the spectral problem for linear acenes [28] . Under re- 
placement £j =>■ k y , Eq. (All) converts into the disper- 
sion relation for polyacene [24]. 

Formally the same equation as Eq. (A10) appears 
in the theory of M-oligomers [60, 61] which are linear 
molecules consisting of N monomers M coupled to each 
other via left and right binding atoms, as illustrated in 
Fig. 8. Monomer M can be described by the Green's 
function G™ a , = (a\(I - H^-^a'), where H M is one- 
particle Hamiltonian in the tight-binding representation. 
Otherwise, M is an arbitrary complex of Nm atoms. 

In Eq. (A10), g 3 a a , associates with the Green's func- 
tion matrix element of a hypothetical monomer indicated 
by the dashed frame in Fig. 8. Thus, all relations which 
follow from Eq. (A10), can be exploited in the present 
context. In general terms, the description of M-oligomers 
with the use of Eq. (A10) represents a generalization of 
the Lennard-Jones theory of polyenes, (M)jv, Nm = 2, 
M = C=C [21]. 

The 7r electron spectrum of conjugated M-oligomers 
is determined by two equations [60, 61]. One of them 
relates the state energy E with the wave vector expressed 
in units of periodicity of the oligomer chain, 



f(G M ) 



(A12) 



where f(G M ) is a functional of monomer Green's function 
matrix elements, 



/(G) = 



(l + Gf r 2 



(A13) 



In our case, k = V3k x , k x is in units of a 1 . The other 
equation determines allowed values of the wave vector 
which must satisfy equation 



sin kN 



G 



M 

l,r 



sian(N + l) Gf t 



M2 



r<M2- 



(A14) 



Equations (4) and (5) follow from Eqs. (A13)-(A15) after 
replacing Gf\ and Gf^. by g\ l and gj r and some algebra. 

It is easy to see that the analogy with M-oligomers is 
straightforwardly applicable to segments of zigzag carbon 
tubes. The only difference is that in Eqs. (A10) and 
(All), we have to replace = jft. by = ^f-, j = 
0, 1, ...,J\f—l. This is to say that the role of hypothetical 
monomer (framed box in Fig. 8) is played by cyclacene 
chain (instead of linear acene chain). 

In this case, applying Eqs. (A12)-(A14), one obtains 



E± 2 = 1 ± 4 



cos 



7TJ 



cos 



4 cos 



2 

AT 



(A15) 



where for the sign + and — , and each value of j, except 
j = N/2, N values of k- , n- v ., Vj = 0,1, ...,N-l, are 
solutions to 



sin nfN 



sin nf(N + 1/2) 



T2 



cos 



7TJ 



for j < [AT/2], and 



sin/tj N 



sin K f(Ar + l/2) 



±2 



cos 



7TJ 



(A16) 



(A17) 



for j > [Af/2]. Here, [x] is the nearest integer function. 

Equations (4) and (A15) were used in Sec. 2 to obtain 
dispersion in graphene daughter structures. The wave 
functions for these structures are discussed next. 



APPENDIX B: WAVE FUNCTIONS [62] 

The wave function (2), which corresponds to the energy 
Ej(Kj), can be represented as 



E 

n a=l,r,X,p 



\m, n, a), (Bl) 



where A - K ± is the normalization constant, 



l 3,K i — / 



^^(Ar + l-n)], 



-s\ sin 
sin 

± sin2teipsm[«±(n-l/2)], 

7Tj(m-l/2) „;„ r„± I 



T*i sin sm[Kf(N + l/2-n)], 

(B2) 

Si = sign [siriKj(iV" + 1/2)], and a — l,r,X,p from up to 
down. For each j, there is N solutions of Eq. (5) with 
plus and minus in its right hand side. 

Thus defined, coefficients (B2) satisfy boundary con- 
ditions (A5). For segments of armchair carbon tubes 
(aCTS), we have, instead of Eq. (A5), 



ipo,; 



IpO,; 



1pm,l,a — 1pm,N+l,a- 



(B3) 
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With these boundary conditions, it is convenient to rep- 
resent energies of tt electron states in aCTS as follows 



The vth root is bounded from below by . Therefore, 
we can redefine k" as follows, 



aCTS 

E 3i' 



Z.; 



jj' I > 



(B4) 



TT ±S2e ^™' (B5) 
j = l,2,...,Af,j' =0,1,...,N-1, 



2cOS 2(Af+l) 



and S2 = sign (cos ^-). 

In these notations, eigen states of aCTS take the form 



m n a=l,r,\,p 



(B6) 



where 



*r m,n,< 



±S 2 - 



Z ML sin ^ m 

sin ^ e 2 -» 
(AT+1) 6 

± S2 sin ^(^-l/2) e2t7r £i^) 



sin 



(B7) 

with the same correspondence between a and l,r, X, and 
p as in Eq. (B2). 

The wave functions for other structures discussed in 
the main text can be obtained via obvious modifications 
of Eqs. (Bl) and (B2) or Eqs. (B6) and (B7). 



APPENDIX C: PROOF OF EQUATION (30) 



V3A£ = -£(i/ + A„). (C3) 

It can be proved that < A„ < 1. Moreover, it turns 
out that for sufficiently small v « N, A„ w |. In this 
approximation, the substitution of (C3) into the left hand 
side of Eq. (20) leads to 



sm 



sin V3k»{N+ 1/2) 



ir(v + A„) 



2N 



cotan(7rA„) 



On the other hand, for q « 1, 



(C4) 



2cos^ w 1- V3q/2-q 2 /8. (C5) 



By equating these two expressions of the same quantity, 
we obtain 



tan(7rA„) 



tt(v + A„) 
V3giV 



(C6) 



The latter equation can be solved approximately, thus, 
solving the initial problem. 

tt(^+1/2) 



Let us assume that \q\ < ■ Then, 



n{v + l/2) 



V3iV n{v + 1/2) ' 



(C7) 



Equation (20) represents the relation between 
cos(fc y /2) and k x in the graphene sheet framed by two 
armchair-shaped and two zigzag-shaped boundaries, as 
shown in Fig. 1. The substitution of this specific expres- 
sion of cos(fcj,/2) into Eq. (19) and usage of a chain of 
trigonometric identities yields 



- i/ y/3kx 

1 - 4 cos — cos 

2 2 

sin 

sin V3k x (N+ 1/2)' 



4cos 2 ^ 
2 



(CI) 



This gives (24). 

All N roots of Eq. (20), < k v x < tt, v = 0, 1, . . . , N - 
1, satisfy the following inequality 



sinV3fc£iV 
< =— — - < 2. 

smV^K{N + 1/2) 



(C2) 



showing that real solutions of Eq. (20) depend linearly 
on the wave vector q « 1. This is in contrast with 
the exponential dependence of the imaginary solutions 
discussed in Sec. HID. If q < q c , 



tt(i/ + 1/2) 



(C8) 



Equations (C7) and (C8) can also be obtained from 
Eq. (17). 

These results prove Eq. (30) and justify a simple 
picture of zGR spectrum above and below the band of 
edge states (but not far from the Fermi energy), where 
q and k x can be considered as independent quantum 
numbers, particularly, near the bottoms/tops of conduc- 
tion/valence bands. 



[1] P. R. Wallace, Phys. Rev. 71, 622 (1947). 



[2] S. Iijima, Nature (London) 354, 56 (1991). 



13 



S. Iijima, T. Ichihashi, and Y. Ando, Nature (London) 
356, 776 (1992). [31 
S. Iijima and T. Ichihashi, Nature (London) 363, 603 [32 
(1993). 

D.S. Bethune, C.H. Kiang, M.S. de Vries, G. Gorman, [33 
R. Savoy, J. Vazquez, and R. Beyers, Nature (London) 
363, 605 (1993). [34 
K.S. Novoselov, A. K. Geim, S.V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, [35 
Science 306, 366 (2004). [36 
K.S. Novoselov, D. Jiang, T. Booth, V.V. Khotkevich, 
S.V. Morozov, and A. K. Geim, Proc. Nat. Acad. Sc. [37 
102, 10451 (2005). 

K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, [38 
M.I. Katsnelson, I.V. Grigorieva, S.V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). [39 
Y. Zhang, Y. -W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). [40 
A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). [41 
R. Saito, M. Fujita, G. Dresselhaus, and M. S. Dressel- [42 
haus, Phys. Rev. B 46, 1804 (1992). 

J.W. Mintmire, B.I. Dunlap, and C. T. White, Phys. [43 
Rev. Lett. 68, 631 (1992). 
K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dres- [44 
selhaus, Phys. Rev. B 54, 17954 (1996). 

R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Phys- [45 
ical Properties of Carbon Nanotubes (imperial College, [46 
London, 1998). " [47 
J.W. Mintmire and C.T. White, J. Phys. Chem. B 109, 
52 (2005). [48 
Y. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Int. 
J. Mod. Phys. B 21, 4611 (2007). 

K. N. Houk, P. S. Lee, and M. Nendel, J. Org. Chem. 66, [49 
5517 (2001). 

R. Pariser, J. Chem. Phys. 24, 250 (1956). [50 

R. Pauncz, J. de Heer, and P.O. Lowdin, J. Chem. Phys. 

36, 2247 (1962). [51 

K. Tanaka, s. Yamashita, H. Yamabe, and T. Yamabe, 

Synthetic metals 17, 143 (1987). [52 

J. E. Lennard- Jones, Proc. Roy. Soc. Ser. A 158, 280 [53 

(1937). 

H. C. Longuet-Higgins and L. Salem, Proc. Roy. Soc. Ser. [54 
A 251, 172 (1959). 

W.P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. B [55 
22, 2099 (1980). [56 
S. Kivelson and O.L. Chapman, Phys. Rev. B 28, 7236 [57 
(1983). 

L. Malysheva and A. Onipko, Phys. Rev. Lett. 100, [58 
186806 (2008). 

L. Malysheva and A. Onipko, arXiv:0802.1385v2 [cond- [59 
mat.mes-hall]. 

L. Malysheva and A. Onipko, arXiv:0803.1761vl [cond- [60 
mat.mes-hall]; JETP, 134, issue 6(12) to appear. 
L. Malysheva and A. Onipko, phys. stat. sol. (b), 245, [61 
2132 (2008). 

Yu. Klymenko, L. Malysheva, and A. Onipko, phys. stat. 
sol. (b), 245, 2181 (2008). [62 
S. Reich, C. Thomsen, and J. Maultzsch Carbon Nan- 
otubes: Basic concepts and Physical Properties (WILEY- 



VCH, Weinheim, 2004). 

T. Ando, J. Phys. Soc. Japan 74, 777 (2005). 

J.-C. Charlier, X. Blase, and S. Roche, Rev. Mod. Phys. 

79, 677 (2007). 

H.-Y. Zhu, D.J. Klein, T. G. Schmalz, A. Rubio, and 
N.H. March, J. Phys. Chem. Solids 59, 417 (1998). 
J.W. Mintmire and C.T. White, Phys. Rev. Lett. 81, 
2506 (1998). 

K. Kobayashi, Phys. Rev. B 48, 1757 (1993). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusak- 

abe, J. Phys. Soc. Japan 65, 1920 (1996). 

K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, 

Phys. Rev. B 59, 8271 (1999). 

Y.-W. Son, M.L. Cohen, and S.G. Louie, Phys. Rev. 
Lett. 97, 216803 (2006). 

C. T. White, J. Li, D. Gunlucke, and J.W. Mintmire, 
Nano Letters 7, 825 (2007). 

D. Gunlycke and C. T. White, Phys. Rev B 77, 115116 
(2008). 

L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 
H. Zheng, Z. F. Wang, T. Luo, Q. W. Shi, and J. Chen, 
Phys. Rev. B 75, 165414 (2007). 

M. Kohmoto and Y. Hasegava, Phys. Rev. B 76, 205402 
(2007). 

A.M. Akhmerov and C.W.J. Beenakker, Phys. Rev. B 
77, 085423 (2008). 

Proof of Eqs. (30) and (31) belongs to Lyuba Malysheva. 
L. Malysheva and A. Onipko, unpublished. 
A. Messiah, Quantum Mechanics (NHPC, Amsterdam, 
1969), Vol 2. 

Dispersion of edge states within the energy interval \E\ < 
(2N + 1) _1 is not covered by Eq. (32) and has to be 
considered separately as explained in Sec. 3D. 
T. Ando and T. Nakanishi, J. Phys. Soc. Japan 67, 1704 

(1998) . 

V.V. Cheianov and V.I. Fal'ko, Phys. Rev. B 74, 
041403(R) (2006). 

M. I Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 

Physics 2, 620 (2006). 

O. Klein, Z. Phys. 53, 157 (1929). 

N. Dombey and A. Calogeracos, Phys. Rep. 315, 41 

(1999) . 

R. Landauer, IBM J. Res. Dev. 1, 233 (1957); 32 306 
(1988). 

M. Biittiker, Phys. Rev. Lett. 57 1761 1986. 
M. Biittiker, Phys. Rev. 38 9375 1988. 
Reference [56] contains a compressed review of the ap- 
proach and lists related publications. 
N. M. R. Peres, A. H. Castro Neto, and F. Guinea, Phys. 
Rev. B 73, 195411 (2006). 

C. W. J. Beenakker and H. van Houten, Solid State Phys. 
44, 1 (1991). 

A. Onipko, Yu. Klymenko, and L. Malysheva, J. Chem. 
Phys. 107, 5032 (1997). 

A. Onipko and L. Malysheva, Nano and Molecular Elec- 
tronics Handbook (CRC Press, Boca Raton FL., 2007), 
Chapter 23. 

Appendixes B and C are authored by L. Malysheva (un- 
published) . 



